function s = getAsymmetricGaussianScore(score_, mean_, sigma_lo_, sigma_hi_)

   

    sigma = 1.0;

    if (score_ == mean_)
        s = 1.0;
    elseif (score_ > (mean_ + 3*sigma_hi_))
        s = 0.0;
        %disp('too hi')
    elseif (score_ < (mean_ - 3*sigma_lo_))
        s = 0.0;
        %disp('too lo')
    else
        if (score_ > mean_)
            sigma = abs(sigma_hi_);
        elseif (score_ < mean_)
            sigma = abs(sigma_lo_);
        end
        
        z = (score_ - mean_) / sigma;
        s = exp(-(z^2.0)/2.0);
        
    end
    
     %a = [score_ mean_ sigma_lo_ sigma_hi_ s];
     %a
    
end